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We studied single W boson production in high energy ej collisions and the sen- 
sitivity of various observables to the WWj gauge boson coupling by evaluating the 



Qh! helicity amplitudes of all Feynman diagrams which contribute to the final state being 

studied, including the W decay to final state fermions. We examined W production 



at 500 GeV and 1 TeV e + e colliders, comparing results for photon spectra obtained 
from a backscattered laser and from beamstrahlung radiation. Here we found that 
the couplings could best be measured using the backscattered laser photons with 
\Sk 7 \ < 0.07 and |A 7 | < 0.05 at a 500 GeV collider and |Jk 7 | < 0.07 and |A 7 | < 0.02 
at a 1 TeV collider, all at 95% C.L.. The measurement of k 7 is at the threshold of 
being able to measure loop contributions to the trilinear gauge boson vertex. For 
completeness we include the limits achievable using single W production at a 200 GeV 
e + e~ collider in the Weizacker- Williams approximation. This process can measure 
K-y to ±0.15 at 95% C.L. which is comparable to the W pair production process. 
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I. INTRODUCTION 



Despite the fact that the standard model of the electroweak interactions [|T| agrees ex- 
traordinarily well with all existing measurements [fJH there is a widespread conviction that 
it is nothing more than a low energy limit of a more fundamental theory ||. An approach 
which is receiving considerable attention is to represent new physics by additional terms 
in an effective Lagrangian expansion and then to constrain the coefficients of the effective 
Lagrangian by precision experimental measurements ||,|5]Hn| • The bounds obtained on the 
coefficients can then be related to possible theories of new physics. For example, this ap- 
proach has been used by a number of authors to bound dimension four operators which can 
contribute to the vacuum polarization tensors of the massive gauge bosons via a global anal- 
ysis of neutral current data ||. These bounds have put severe constraints on technicolour 
theories of dynamical symmetry breaking. Similarly, the trilinear gauge boson couplings 
have been described by effective Lagrangians ||,[K|[II] . In one commonly used parametriza- 
tion, for on shell photons, the CP and P conserving ^WW vertex is parametrized in terms 
of two parameters, k 7 and A 7 (1(J. Although bounds can be extracted from high precision 
low energy measurements and measurements at the Z° pole [0, there are ambiguities and 
model dependencies in the results [|13j. In contrast, gauge boson production at colliders 
can measure the gauge boson couplings directly and unambiguously. The best current di- 
rect measurement of these parameters come from the CDF and DO collaborations at the 
Fermilab pp. Their combined limits, varying one parameter at a time, are \8kj\ < 0.9, 
|A 7 | < 0.3, \5kz\ < 0.5, and \Xz\ < 0.3 at 95% C.L. |yj. Putting tight constraints on the 
trilinear gauge boson couplings is one of the primary motivations for the LEP200 upgrade 
l^ , p!5lJl6| , p!7|JT8| , ^9|] . At the time of writing the 95 % C.L. bounds from the LEP experiments 



17j are —0.9 ~ 5k ~ 1.12 and —0.78 ~ A ~ 1.19 where only one parameter is varied at 



a time and these results have assumed that 5kj = 5k,z and A 7 = Xz- The L3 collabo- 
ration at CERN has also obtained the 95% C.L. bounds from single W production |18 



-3.6 < 5k 7 < 1.5 and —3.6 < A < 3.6. Precise direct measurement to the level of several 
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percent will have to wait until the era of the Large Hadron Collider |T9|j20[| and the NLC, a 
high energy e + e~ collider. 

Recently, the idea of constructing ej and 77 colliders using either high energy photons 



from lasers backscattered off of a high energy electron beam |HJ or photons arising from 
beamstrahlung radiation |2"2] , |2"3"| , |2~3|] has received serious attention. The physics possibilities 
for e7 colliders are the subject of a growing literature [f25|| . In particular, the properties of W 
bosons, including the ^WW coupling, has been examined in a number of recent publications 
2|.|27|. The LEP collaborations p| are starting to obtain results for single W production 



where the photon arises in t-channel exchange |28[ although the results are still rather weak 
compared to those obtained by CDF and DO. 

In this paper we reexamine the sensitivity to which the 'jWW vertex can be measured 
at e7 colliders using photon spectra produced from backscattered lasers and beamstrahling 
radiation. In contrast to other analyses p7| we include the decay of the W boson to final state 
fermions along with contributions to the final state that do not proceed via an intermediate 
W boson. At high energy this provides additional information off the W resonance through 
the interference of the various diagrams. In addition, we considered the various backgrounds 
that may obscure results using hadronic W decay. 

We focus on W production at 500 GeV and 1 TeV e + e~ colliders and compare the sensi- 
tivities achievable using a backscattered laser photon spectrum and a beamstrahlung photon 
spectrum. For completeness we also include results obtained from single W production at 
a 200 GeV e + e~ collider using the Weizacker- Williams effective photon approximation as a 
competing process to W pair production. 

In the next section we write down the effective vertex and the resulting Feynman rule. In 
Sec. [TTT] we present our calculation and results and in the final section we give our conclusions. 
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II. THE W W7 EFFECTIVE VERTEX 



Within the standard model the WW'j vertex is uniquely determined by SU(2)l x U(l) 
gauge invariance so that a precise measurement of the vertex poses a severe test of the 
gauge structure of the theory. The most general WW'j vertex, consistent with Lorentz 
invariance, can be parametrized in terms of seven form factors when the W bosons couple 



to essentially massless fermions which effectively results in d^W^ = |ID|,n|. For on shell 
photons, electromagnetic gauge invariance further restricts the tensor structure of the WW^f 
vertex to allow only four free parameters, two of which violate CP invariance. Measurement 
of the neutron electric dipole moment constrains the two CP violating parameters, k 7 and 
A 7 , to values too small to give rise to observable effects in the process we are considering; 



k~\, |AJ < O(10 4 ) |29, 30|. Therefore, the most general Lorentz and CP invariant vertex 



compatible with electromagnetic gauge invariance is commonly parametrized as | lQfl : 



C WWl = -le I (Wj^M* - WlA u wn + k^W^ + ^W^F**} (1) 

where A 11 and W^ are the photon and W~ fields, = d^W u -d u W^ and F^ v = d^Ay-dyA^ 
denote the W and photon field strength tensors, and Mw is the W boson mass. Higher 
dimension operators correspond to momentum dependence in the form factors. The first 
term, referred to as the minimal coupling term, simply reflects the charge of the W. The 
Feynman rule for the WW'j vertex resulting from eqn. (1) is given by: 

ie{ g a p[(l — Xk-- q)k +fl - (1 - A k+ ■ q)k- fl ] 

- g af i[(l -Xk-- q)k +fj - (k - A fc+ ■ fc_)g (U ] 

- 9pn[( K - Xk ~ ' k +)^ ~(l-Xk + - q)k_ a ] 

+ X(k +fl k- a q/3-k-^q a k +/3 )} (2) 

with the notation and conventions given in Fig. 1 and where A = A/Mjy. At tree level the 
standard model predicts k 7 = 1 and A 7 = 0. Other parametrizations exist in the literature 
such as the chiral Lagrangian expansion and one can map the parameters we use to those 
used in other approaches ||[|. 



Figure 1: Feynman diagram for the WWj vertex corresponding to the Lagrangian and Feynman 
rule given in the text. 

Before proceeding it is useful to comment on what to expect for k 7 and A 7 . From the 
chiral Lagrangian approach one expects that <5k 7 ~ O(10~ 2 ) and A 7 is suppressed by an 
additional factor of 100 ||. These order of magnitude estimates are confirmed by explicit 
calculation. The contribution of the t-quark results in 5k ~ 1.5 x 10~ 2 and Ay — +2.5 x 10 -3 



while a SM Higgs boson of mass 200 GeV contributes 5k ~ 5 x 10 and A ~ 4 x 10 [31 



Technicolour theories give 5k z = —0.023 and 5k^ = [0 and supersymmetric theories give 
SK-max ^ 7 x 10 -3 and X max ~ 10 -3 |[J2|. 

What one gleans from these numbers is that a deviation of more than a couple of percent 
would be difficult to accomodate in the SM and in general, contributions via loop corrections 
typically contribute no more than a couple of percent. Although it might be possible to find 



models that give slightly larger contributions, for example models with Z — Z' mixing [33 



or models with many particles which contribute coherently to loop contributions |34j], if 
deviations much larger than several percent are observed this would signal something very 
radical such as composite gauge bosons ||35|| . Since we know of no convincing models of this 
sort, to probe for new physics via anomalous trilinear gauge boson couplings one must be 
able to measure the vertex to the level of a few percent. 

Deviations from the standard model (a = 5k — k — 1, A) lead to amplitudes which grow 
with energy and therefore violate unitarity at high energy |36| , |3~T|| . One method of avoiding 
violation of the unitarity bound is to include a momentum dependence in the form factors, 
a (<?w/; ^7 = 0)' so that the deviations vanish when either \q^\ or the absolute square 



of the four momentum of the vector bosons, becomes large [p7|| . We therefore include the 
form factors 
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a(q 2 w , f w , 0) = ao[(l + |^|/A 2 )(1 + \f w \/A 2 )}- 



n 



(3) 



where A represents the scale at which new physics becomes important and n is chosen as the 
minimum value compatible with unitarity. We take n — 1 and A = 1 TeV in our numerical 
results. We find that lowering this scale effects our conclusions slightly although in this case 
the new physics should show up elsewhere. Increasing the scale has only a small effect on 
our results. 



The WW'y vertex we are studying contributes via diagram 2b. To preserve electromagnetic 
gauge invariance and to properly take into account the background processes our calculation 
includes all the diagrams of Fig 2 for arbitrary values of k 7 and A 7 . To obtain the cross 
sections and distributions we used the CALKUL helicity amplitude technique which for 
completeness we summarize in Appendix A along with the amplitudes corresponding to Fig. 
2. Monte Carlo integration techniques are then used to perform the phase space integrals 
f39|j . We treat the photon distributions as structure functions, / 7 / e (a;) and integrate them 
with the cj cross sections to obtain our results: 



where x is the fraction of the original electron energy carried by the photon. For completeness 
and for the convenience of the interested reader we include the various photon distributions 
in Appendix B. For our numerical results we take a(M z ) = 1/128, M w = 80.22 GeV, 
T w = 2.0 GeV, sin 2 6 W = 0.23. 



III. CALCULATION AND RESULTS 



The Feynman diagrams contributing to the process e 7 



vff are given in Fig. 2. 




(4) 
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Figure 2: The Feynman diagrams contributing to the process e^f — ► vqq. For the process — > 
v e ^ii diagram (c) does not contribute and the quark charge in diagram (d) should be replaced by 
the \x charge. 

The signal we are studying consists of either, (i) for leptonic W decay, a high transverse 
momentum lepton (p?) and large missing transverse momentum {£>t) due to the neutrinos 
from the initial electron beam and from the W decay, or (ii) for hadronic W decay, two 
hadronic jets and large missing transverse momentum (for) due to the neutrino from the 
initial electron beam. In both cases, we impose the acceptance cut that visible particles in 
the final state be at least 10° from the beam direction. Our conclusions are not sensitive to 
the exact value of this cut. We also impose a cut on the minimum £>t] for y/s = 200 GeV 
we used /p T > 5 GeV and for ^/Ts = 500 GeV and 1 TeV we used /p T > 10 GeV. We do 
not include fragmentation and hadronization effects for the hadronic modes and identify the 
hadron jet momenta with that of the quarks. The signals we consider are therefore 

e~+7^^-+j6 (5) 
e~ +7 ->3+3+ j> (6) 

We also examined the reaction e~7 — * evv which includes the process e~7 — > e~Z° — > e~vv 
in addition to the diagrams of figure 2. However, once kinematic cuts are imposed that 
eliminate the uninteresting contributions from — > e~Z we are left with results comparable 
to the n mode. 
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A. Backgrounds 



Before proceeding we discuss potential backgrounds to the processes we are studying 
PUp . They can be divided into 3 categories: 

1. Direct 

e~7 — > e~ Z° — > e~gg (7) 
e~e + — > 7Z — > 755 (8) 

where in the first case the outgoing e~ is not observed while in the second case the 
outgoing 7 is not observed. 

2. Once Resolved. These involve 79 or 7$ collisions where the q or g is treated as a parton 
in the photon. The subprocesses are: a{piq — > qg), cr(-yq — > gg), and (7(7(7 — > qq). 

3. Doubly Resolved. Here, q or g partons come from both the 7 and a Weizacker- Williams 
photon coming from the electron. The subprocesses are: a(gg — > qq), &{gg — > gg), 

°(Q9 ~> <%9 00), %9 ??), 5V), <%9 99), &(q<f -> 9<f)> and 

a(gg' -> gg'). 

The e~e + — > Z°7 is easily removed by imposing the constraint that the photon not be 
observed while at the same time there is significant missing pt in the event. Since the photon 
recoils against the Z these two contradictory requirements eliminate this background. We 
show the jet — jet invariant mass for the signal and remaining backgrounds for the backscat- 
tered laser case with y/s = 500 GeV in Fig. 3. These results include the detector acceptance 
cuts of 10° < Q e -jet < 170° and pr(jei) > 5 GeV. The resolved photon backgrounds are 
substantial but they can be eliminated by requiring a minimum missing pj>. This is a simple 
consequence of the fact that we require both final state particles to be observed while the 
initial state partons are essentially collinear with the beam direction. The only remaining 
background is the direct background e~7 — ► eZ with the electron escaping down the beam 
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pipe. We eliminate this background by imposing a tight cut on the invariant mass of the 
outgoing jets so only jets contributing to the W peak remain. This eliminates the bulk of 
the background from the large Z peak. For the case of dijet pairs with invariant mass off 
the W resonance we can eliminate the e~7 — > eZ background by imposing a much larger 
cut on the missing p T . Since we eliminate events when the electron is observed no events 
will remain with missing p F greater than the maximum pt c for an electron lost down the 
beam pipe (~ 40 GeV for y/7s = 500 GeV). 
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Figure 3: Dijet invariant mass for the signal and backgrounds. The solid line is the signal, while 
the dashed line is the signal with the detector resolution p T cut. The short dashed curve is the 
e7 — > eZ background, the long dashed curve is the total of the singly resolved backgrounds and 
the dotdashed curve is the total of the doubly resolved backgrounds. 



B. Results 

In the following subsections we study single W production for y/s = 200 GeV, 500 GeV, 
and 1 TeV. In doing so we examined a number of kinematic distributions, asymmetries, and 
ratios of observables for their sensitivity to anomalous gauge boson couplings. Since we will 
only present results for the observables most sensitive to the gauge boson couplings, at this 
point we list all observables that we studied. For the leptonic W decay modes we considered; 
da/dpT^L, da/dE^, da/dcos9 efi , where 9 e)M is the angle between the incoming electron and the 
outgoing muon, A FB = (o> — <J B )/(cr F + a B ), R IO = <Ti/<t where 07 is the cross section for 
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I cos# eM | < 0.4 and a is the total cross section with the kinematic cuts given above. When 
backscattered laser photons are used we can form the ratio R13 = <J\ja^ where <j\ is the 
cross section for the mainly helicity 1/2 amplitudes which result when the incident photons 
are left handed and 03 is the cross section for the mainly helicity 3/2 amplitudes when the 
incident photons are right handed^. For the hadronic W decay modes we reconstructed 
the W boson 4-momentum from the hadronic jets' 4-momentum, imposing the kinematic 
cut of 75 GeV < M qq = \J (p q + p q ) 2 < 85 GeV. Including the nonresonant diagrams of fig 
(2c) and (2d) and reconstructing the W boson in this manner gives different results than 
from simply studying the cross sections to W bosons. Finally, in some cases, which we will 
describe in further detail below, we studied the effect of anomalous couplings on the cross 
section of the hadronic modes off the W pole (i.e. M qq 7^ Mw). In general, deviations 
of the gauge boson couplings had a very substantial effect on the cross section off the W 
resonance although, because of the reduced cross section, the statistical significance is not 
really enhanced. This does, however, point out the importance of considering the process 
that is actually measured, not just a theorist's idealization. 



C. yfs = 200 GeV: LEP2 

This process has been examined in detail by a number of authors and the various LEP 
collaborations have presented preliminary results. As our results are consistent with other 
work, we include the limits that are achievable at LEP2 mainly for completeness. At LEP 
200 the possibility of cy collisions is only possible via t-channel photon exchange between 
the electron and positron beams | 28f . We therefore examine single W production in e + e~ 



collisions by folding the cross sections for ey — > qqv and ey — > fivv with the Weizacker- 
Williams photon distribution. 

We find that distributions utilizing W's reconstructed from the hadronic modes are most 



1 Note that the high energy photon beam has opposite polarization to that of the laser. 
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sensitive to the WW^y couplings. This is primarily due to the increased statistics of the 

hadronic modes over the leptonic modes since the bounds obtained at 200 GeV are limited 

by statistics; = 0.059 pb vs aw = 0.28 pb (with the cut on M q q). To obtain sensitivity 

bounds we used the total cross section for the reconstructed W boson. The 95% C.L. limits 

are given in Fig. 4. To obtain this curve we calculated the observable for a given value of k 7 

and A 7 and determined at what level the value, if measured, would be compatible with the 

standard model prediction. The 95% C.L. limits obtained by varying one parameter at a 

time, for L = 500 pb" 1 , are 0.85 < k 7 < 1.13 and -0.63 < A 7 < 0.61 from a(e + e~ -> e + W~). 

The sensitivity on k 7 is comparable to that obtained from the TV-pair production process 

but because it offers a means of measuring the WW'y vertex independently of the WW Z 

vertex, it is an important complement to the PF-pair production process. 
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Figure 4: The achievable bounds on k 7 and A 7 at 95% C.L. for a y/s = 200 GeV e + e~ collider using 
the total cross section for the hadronic decay mode of the W boson. Fig. (a) is for L = 100 pb -1 
and fig. (b) is for L = 500 pb -1 . 



D. The Next Linear Collider; </s = 500 GeV and ^ = 1000 GeV 



For the NLC [41] we consider two possibilities for the photon spectrum; that arising from 
backscattering a laser from one of the original electron beams and beamstralung, which is the 
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radiation which arises when two intense beams of electrons pass through one another. For 
the results using the beamstrahlung photon spectrum we concentrate on the beam spectrum 
resulting from the G set of parameters of Ref. [Q. We will discuss the effect of different 
beam parameters on our results. We included the Weizacker- Williams contributions in our 
beamstrahlung results. 

The NLC is envisaged as a very high luminosity collider so that the number of events 
per unit of R, the QED point cross section, which is an s-channel process and goes like 
1/s, remains reasonable. The integrated luminosities for a Snowmass year (10 7 sec) are 
expected to be ~ 60 fb _1 for a \/s = 500 GeV collider and ~ 200 fb _1 for a y/s = 1 TeV 
collider. Typical cross sections for the process — > fiv^Ve and ey — > Wv — > qqv at 
\fs = 500 GeV are 3.2 pb and 16.6 pb respectively for the backscattered laser mode and 
2.1 pb and 10.8 pb respectively for the beamstrahlung mode, leading to ~ 10 6 events per 
year. At \/s = 1 TeV a(e r y — > fiv^v e ) = 4.0 pb and a(e r y — > Wv — > qqv) = 19 pb for the 
backscattered laser mode and a^e-y — > fiv^v e ) = 6.0 pb and a(ej — > Wv — > qqv) = 31 pb 
for the beamstrahlung mode leading to ~ 6 x 10 6 events/year. Thus, except for certain 
regions of phase space, the errors are not limited by statistics, but rather by systematic 
errors. Factors entering into systematic errors include an accurate knowledge of the total 
luminosity, particle misidentification, triggering and detector efficiencies, uncertainty in the 
size of backgrounds, calorimetric accuracy etc. Estimating systematic errors requires detailed 
Monte Carlo studies which we do not attempt. For cross sections we assume a systematic 
error of 5% and for asymmetries and ratios, where some of the systematic errors cancel, we 



assume a systematic error of 3% j43fl . We consider the effects of reducing these systematic 
errors on the achievable sensitivities. In our results we combine in quadrature the statistical 
errors based on the integrated luminosities given above with the systematic errors: 

52 = S stat + 5 sys- ( 9 ) 

We begin with y/s = 500 GeV. The total cross sections and the angular distributions of 
the outgoing muon and reconstructed W (from the hadronic decay mode) are sensitive to 
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anomalous couplings. We plot the distributions for both the backscattered photon case and 
the beamstrahlung case in Fig. 5. At higher energies we can obtain additional information, 
especially for A 7 , from the spectrum of the outgoing lepton or the reconstructed W. We 
show these spectra in Fig. 6. Finally, as already pointed out, the invariant mass distribution, 
shown in Fig. 7, of the qq pair above the W mass also provides useful information. If, for 
example, we integrate the M qq - spectrum from 100 GeV up, we obtain a cross section of 0.25 
pb for the backscattered laser mode which offers considerable statistics. For M qq > 300 GeV, 
a = 0.006 pb which yields ~ 400 events/year. More importantly, this high M qq region shows 
a higher sensitivity to anomalous couplings than the M\y pole region. 
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Figure 5: The angular distributions of the outgoing muon and reconstructed W boson relative 
to the incoming electron for *Js = 500 GeV. (a) For a muon with the backscattered laser photon 
spectrum, (b) for a muon with the beamstrahlung photon spectrum, (c) for a reconstructed W 
boson with the backscattered laser photon spectrum, and (d) for a W boson with the beamstrahlung 
photon spectrum. In all cases the solid line is the standard model prediction, the long-dashed line 
is for k 7 = 0.6, A 7 = 0, the short-dashed line is for re 7 = 1.4, A 7 = 0, and the dotted line is for 

K-y = 1, A 7 = 0.4. 
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Figure 6: The pt distributions of the outgoing muon and reconstructed W boson for yfs = 
500 GeV. (a) For a muon with the backscattered laser photon spectrum, (b) for a muon with 
the beamstrahlung photon spectrum, (c) for a reconstructed W boson with the backscattered laser 
photon spectrum, and (d) for a W boson with the beamstrahlung photon spectrum. In all cases the 
solid line is the standard model prediction, the dotted is for k 7 = 0.6 and A 7 = 0, the long-dashed 
line is for k 7 = 1.4, and A 7 = 0, and the dot-dot-dashed line is for k 7 = 1, A 7 = 0.4. 
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Figure 7: The hadron jet invariant mass (M q q) distribution for y/s = 500 GeV. (a) For the 
backscattered laser photon spectrum and (b) for the beamstrahlung photon spectrum. In both 
cases the solid line is the standard model prediction, the long-dashed line is for k 7 = 0.6, A 7 = 0, 
the short-dashed line is for k 7 = 1.4, A 7 = 0, and the dotted line is for k 7 = 1, A 7 = 0.4. 

To quantify the observations of the above paragraph we consider the cross sections, the 
angular distributions, and the distributions for the muon and hadronic modes. For the 
angular distributions we used four equal bins and for the distributions we used the 4 p T 
bins; - 100 GeV, 100-150 GeV, 150-200 GeV, and 200-250 GeV. The bounds obtained for 
these observables are shown in Fig. 8 based on systematic errors of 5%. Bounds obtained 
by varying one parameter at a time are summarized in Table I. In general the limits on 
/t 7 obtained using the beamstrahlung spectrum are comparable to those obtained from the 
backscattered photon spectrum with the exception that we do not show results using angular 
distributions in the beamstruhlung mode because of the ambiguity in deciding the angle 
due to the poor charge determination of VF's decaying hadronically. If the muon charge 
is identified in the beamstrahlung case, the angular distribution method could be used as 
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described, although we do not show these results in Fig 8 and 11. The bounds on A 7 are 
tighter using the backscattered laser photons which reflects the harder photon spectrum in 
this case to which A 7 is more sensitive. We find that « 7 can be measured to within 7% and A 7 
to within ±0.05 at 95% C.L., using the backscattered laser approach, which is approaching 
the sensitivity required to observe the contributions of new physics at the level of loop 
corrections. One can also obtain limits using the high invariant mass region. A priori one 
might expect single W production to be overwhelmed by the — > Ze background. However, 
as pointed out in a previous section, this background can be eliminated in the kinematic 
region of interest with an appropriate kinematic cut, jpr > 40 GeV, which has little effect 
on the signal. We find that for M q g > 100 GeV we obtain the 95% C.L. sensitivity of 
-0.09 < <5/t 7 < 0.08 (-0.1 < 5k 7 < 0.10) and -0.1 < A 7 < 0.16 (-0.16 < A 7 < 0.24) for the 
back scattered (beamstrahlung) cases. 

1.10 
1.05 
a 1.00 
0.95 
0.90 

-0.2 -0.1 -0.0 0.1 0.2 

X 

1.10 
1.05 
y 1.00 
0.95 
0.90 

-0.2 -0.1 -0.0 0.1 0.2 

X 

Figure 8: The achievable bounds on k 7 and A 7 at 95% C.L. for a y/s = 500 GeV e + e _ collider 
(a) using a backscattered laser photon spectrum (b) using a beamstrahlung photon spectrum. The 
dashed line is based on the angular distribution divided into four bins, the dotted line is based on 
the pr distribution of the W boson divided into the four bins given in the text, and the solid line 
is the combined angular and px bounds. 
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TABLES 

TABLE I. Bounds on k 7 and A 7 from the processes e~j — > yTv^Vf, and e~j — > VL~ — > e + gg 
at = 500 using backscattered laser photon distributions (laser) and beamstrahlung photon 
distributions (beam). 
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-0.20 


+0.21 
-0.21 




beam 


+0.04 
-0.05 


+0.07 
-0.08 


+0.09 
-0.09 




weak bounds 




cos 9w 


laser 


+0.05 
-0.05 


+0.06 
-0.07 


+0.07 
-0.07 


+0.09 
-0.09 


+0.10 
-0.10 


+0.10 
-0.11 




laser 


+0.05 
-0.05 


+0.06 
-0.07 


+0.07 
-0.08 


+0.04 
-0.04 


+0.05 
-0.05 


+0.05 
-0.05 




beam 


+0.05 
-0.06 


+0.06 
-0.07 


+0.07 
-0.08 


+0.05 
-0.05 


+0.05 
-0.05 


+0.06 
-0.06 
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Unlike the case at lower energies, the limits from the muon mode and reconstructed W 
mode are comparable. This is due to two reasons: First, while the W mode is restricted 
to the small portion of the phase space at the W mass. In contrast, the muon mode 
reflects the entire kinematically allowed region, in particular the highest energy region where 
deviations from the standard model are most pronounced. Although the diagram we are 
interested in does not dominate in this higher energy region, the interference between it 
and the non-resonant diagrams are important. Once again, this underlines the importance 
of considering all contributions to the process that will actually be observed and only then 
impose constraints. Furthermore, because of the large expected integrated luminosities the 
errors are dominated by systematic errors so that the differences between the hadronic and 
leptonic branching fractions become unimportant. 

We repeat the above exercise for y^s = 1 TeV using an integrated luminosity of 200 fb _1 . 
As was the case at 500 GeV we find that the angular and px w distributions give the best 
constraints on anomalous couplings. These distributions are shown in Fig. 9 and 10 and are 
seen to be qualitatively similar to the 500 GeV distributions. The invariant mass distribu- 
tions of the qq pair is very similar to the 500 GeV except that it extends out about a factor of 
two further. To extract bounds from these distributions we divided the angular distributions 
into four equal bins and the pr w distribution into the four bins; — 200 GeV, 200-300 GeV, 
300-400 GeV, and 400-500 GeV. At 1 TeV there is little improvement on the sensitivity to 
k 7 and about a factor of two improvement on the sensitivity to A 7 resulting in a possible 
limit of <5A 7 = ±0.016 at 95% C.L. using the laser mode. This measurement improvement 
on A 7 reflects the greater sensitivity at higher energy. Confidence level curves are shown in 
Fig. 11 with some of the results summarized in Table II. For M q g > 600 GeV we obtain 
the 95% C.L. sensitivity of |A 7 | < 0.03 (|A 7 | < 0.05) for the back scattered (beamstrahlung) 
cases. The k 7 sensitivies are relatively weak and are therefore not given. 
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Figure 9: The angular distributions of the outgoing muon and reconstructed W boson relative 
to the incoming electron for y/s = 1 TeV. (a) For a muon with the backscattered laser photon 
spectrum, (b) for a muon with the beamstrahlung photon spectrum, (c) for a reconstructed W 
boson with the backscattered laser photon spectrum, and (d) for a W boson with the beamstrahlung 
photon spectrum. In all cases the solid line is the standard model prediction, the long-dashed line 
is for k 7 = 0.6, A 7 = 0, the short-dashed line is for k 7 = 1.4, A 7 = 0, and the dotted line is for 




100 200 300 400 100 200 300 400 



P T P T 

Figure 10: Thepy distributions of the outgoing muon and reconstructed W boson for y/s = 1 TeV 
with the same labelling as Fig. 9. 
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Figure 11: The achievable bounds on k 7 and A 7 at 95% C.L. for a ^fs = 1 TeV e + e collider 
with the same line definitions as Fig. 8. 



TABLE II. Bounds on k 7 and A 7 from the processes e~7 — >■ ^i"v^u e and e~7 — > W~ —> e 1 qq 
at 1 TeV using backscattered laser photon distributions (laser) and beamstrahlung photon distri- 
butions (beam). 









5k-j 








5A 7 








68% C.L. 


90 % C.L. 


95 % 


C.L. 


68% C.L. 


90 % C.L. 


95 % C.L. 












v^ = 


1 TeV 








laser 


+0.05 
-0.05 


+0.08 
-0.09 


+0.09 
-0.10 




+0.07 
-0.08 


+0.09 
-0.10 


+0.10 
-0.11 




beam 


+0.05 
-0.05 


+0.08 
-0.08 


+0.09 
-0.10 






weak bounds 






laser 


+0.05 
-0.05 


+0.08 
-0.09 


+0.09 
-0.10 




+0.08 
-0.08 


+0.10 
-0.10 


+0.11 
-0.12 




beam 


+0.05 
-0.05 


+0.08 
-0.08 


+0.09 
-0.10 






weak bounds 




cos 6w 


laser 


+0.05 
-0.06 


+0.06 
-0.07 


+0.07 
-0.08 




+0.04 
-0.04 


+0.05 
-0.05 


+0.05 
-0.05 


PT W 


laser 


+0.05 
-0.07 


+0.06 
-0.10 


+0.07 
-0.11 




+0.013 
-0.013 


+0.015 
-0.015 


+0.016 
-0.016 




beam 


+0.05 
-0.07 


+0.06 
-0.10 


+0.07 
-0.12 




+0.015 
-0.015 


+0.017 
-0.017 


+0.018 
-0.018 
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To obtain these sensitivities to anomalous couplings we made a number of assumptions 
on A, the scale of new physics used in the form factors, the beamstrahlung spectrum, and 
the systematic error. We discuss the effects of varying these parameters starting with the 
systematic error which is relevant to the 500 GeV and 1 TeV cases. Reducing the systematic 
error from 5% to 2% reduces the precision on k 7 roughly proportionately with the systematic 
error, i.e. a factor of two reduction in the systematic error will tighten the limits on k 7 by 
roughly a factor of two. In contrast, the attainable constraints on A 7 does not in general 
improve as much, especially for constraints obtained from the pr w distributions which give 
the tightest of all bounds on A 7 . This is more pronounced for the 500 GeV case than the 
TeV case. We can see the reason for this by refering to figures 5, 6, 9, and 10. Varying k 7 
results in an overall shift in the cross section effecting all regions of phase space while in 
contrast, the effect of A 7 grows larger with increasing pr w - The largest effect is at highest 
p T where the cross section, and hence the statistics are lowest. Thus, at least for the 500 
GeV case, statistical errors still play a role and it indicates that at a TeV they could also be 
important if the NLC does not achieve the large integrated luminosities we have assumed. 

We next consider the effect of using the beamstrahlung spectrum arising from the G=l 
beam geometry. The only place this change has any effect on our results is a slight improve- 
ment on the bounds obtained using the pr distributions. This can be attributed to the effect 
the harder photon spectrum has on the px distribution. 

Finally, we consider the effect of varying A, the energy scale used in the form factors to 
preserve unitarity at high energy. We took A = 1 TeV to obtain our results. Changing A 
to 500 GeV and to 2 TeV has no effect whatsoever on our y/s = 200 GeV bounds. For the 
500 GeV case there is a small decrease in the sensitivity (of the order of a few percent) if we 
decrease A to 500 GeV and virtually no change when we increase it to 2 TeV. This is true for 
all measurements except those involving the pr distribution. This is not unexpected. Since 
the form factor was introduced to suppress anomalous couplings at high energy, if the cutoff 
scale is reduced it is doing what is was introduced to do. If this scale is much larger than 
the characteristic energy scale of the process being studied, increasing it further should have 
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no effect. The exception of the px distributions reflects that these particular observables 
are most sensitive to deviations at high energy. We find similar effects for the y/s = 1 TeV 
collider although here the changes are more pronounced. To demonstrate the effect on the 
Pt distribution we plot in Fig. 12 the pr distribution of the outgoing muon for the three 
values of A with A 7 = 0.1 using the backscattered laser spectrum with y/s = 1 TeV. We do 
not consider the choice of this scale an important one. If it were small enough to make a 
difference in these measurements we would expect the new physics to manifest itself in other 
ways, otherwise, the scale would be large enough not to matter. 




Figure 12: The pr distribution for the fi at y/s = 1 TeV and taking A 7 = 0.1 for different 
values of A, the energy scale in the anomalous couplings form factor. The short dashed line is for 
A = 500 GeV, the solid line for A = 1 TeV, and the dashed line for A = 2 TeV. 

IV. CONCLUSIONS 

We examined single W production in e7 collisions for a number of collider energies and 
sources of energetic photons. In our studies we included the W boson decays to final state 
fermions and other processes which contribute to the same final state. At high energy, the off 
resonance results are important since interference effects between these other diagrams and 
the W production diagrams enhance the significance of anomalous couplings, particularly A 7 . 
Although these effects contribute relatively little to the total cross section, their significance 
in constraining the anomalous couplings can be large, especially at high energies and high 
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luminosities where these effects are statistically significant. 

We found that using single W production in e + e~ collisions at y/s = 200 GeV in the 
effective photon approximation k 1 could be measured to ~ ±0.15 and A 7 to ~ ±0.6 at 95% 
C.L.. The former measurement is comparable to what can be achieved in W pair production 
while the latter is not quite as precise. What makes this process interesting is that it offers 
a means of measuring the WW'j couplings independently of the WW Z couplings which is 
far from trivial in the W pair production process. 

For the high energy, high luminosity, NLC e + e~ collider we considered both a backscat- 
tered laser photon spectrum and a beamstrahlung photon spectrum. For yfs = 500 GeV 
5k 7 ~ ±0.07 and 5A 7 ~ ±0.05 and for ^ = 1 TeV, 5^ ~ ±0.07 and 5A 7 ~ ±0.016 using 
the backscattered laser approach. Using the beamstrahlung photon spectrum is only slightly 
less sensitive. The measurement of k 7 is approaching the level of radiative corrections and 
might be sensitive to new physics at the loop level. On the other hand, it is expected that 
the sensitivity to A would have to be at least an order of magnitude more sensitive to be 
interesting. From our analysis it does not appear that there is any overwhelming advantage 
to go to higher energies to study the trilinear gauge boson couplings using ej collisions. 
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APPENDIX A: HELICITY AMPLITUDES 



In this appendix, we summarize the CALKUL spinor technique and give the helicity 



amplitudes for the process — > uff where //' can be fiu^ or qq'. We do not go into 
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any detail and refer the interested reader to the literature and references therein. We limit 
our discussion to massless fermions and massless external gauge bosons, which apply to our 
problem. The propagators for the fermions and gauge bosons have the same form as in the 
usual trace technique so we will not discuss them here. 

The spinor technique results in reducing strings of spinors and gamma matrices to sand- 
wiches of spinors which can be evaluated easily. In doing so, one makes extensive use of the 
right and left projection operators defined by u± = |(1±7 5 ). One defines two four- vectors, 
&o and ki , which obey the following relations: 

ko ■ ko = 0, k\ ■ k\ = —1, ko • k± — 0. 

and the basic spinors: 



u_(k )u_(k ) = u}_fc 



and 



u + (k ) = fau-iko) 

Note that in the massless limit, one can use u and u to describe both particles and antipar- 
ticles, with the spin sum J2x u x(p)ux(p) = £)■ These two spinors are the building blocks for 
any spinor of lightlike momentum p : 

ftu„ x (k ) 
MP) = /0 r 

Two identities are essential for the reduction of the strings; the spin sum given above and 
the Chisholm identity: 

where A is ±1 and represents the helicity state. These two identities allow one to reduce 
strings of spinors and gamma matrices to sandwiches of spinors. Only two of the four 
possible sandwiches are non-zero: 

S(P1,P2) = M+(pi)Mp 2 ) = ~ S (P2,Pl) 
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and 

t(Pl,P2) = «-(Pi)«+(P2) = 8(p 2 ,Pl)*. 

Once the amplitude has been reduced to a series of factors of s(pi,pj) and t(pk,Pi), the 
expressions can be evaluated by computer. A judicious choice of the four-vectors fcg and ki 
simplifies the evaluation of the s and t terms. For our calculation, we used the definition of 
ref. 

# = (pIpIpIp!) 
a# = (1,1,0,0) 
K = (0,0,1,0) 

to obtain 



s(pi,p2) = (pi + m) 1 - w + % V2) 1 

ypx-pf \jP2-pl 

These forms are ideally suited for programming. When dealing with several diagrams, 
one simply evaluates the amplitudes of each diagram as complex numbers and squares the 
sum of the amplitudes in order to obtain the \amplitude\ 2 . 

To include massless gauge bosons one represents, as usual, the gauge boson by its polar- 
ization vector. Following Kleiss and Sterling we use the definition: 

e»{k) = —L=u x (k)^u x (p) 
y4 p ■ k 

where p^ is any lightlike four- vector not collinear to or k§. The choice of p^ acts as a 
choice of gauge and provides a powerful verification of gauge invariance; it can be shown that 
two different choices of p^ will lead to two expressions that will differ by a term proportional 
to the photon momentum. When dotted into the amplitude, this extra term must vanish 
identically because of gauge invariance. Hence, two different choices of p M must give exactly 
the same answer. If they don't, there is a mistake in the amplitude. Generally, we choose 
p^ to be one of the four- vectors of the problem at hand. 

Using this technique we obtain for the helicity amplitudes corresponding to the Feynman 



diagrams of Fig. 2, using the notation M = ieg Mj J4p> 1 ■ k: 
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M LL = B~( r=\ *(P"> 9) S (^' Pc) *(Pe, P 7 ) S ( /C ' Pe) ( Al ) 

s %W + gj 
2 

m lr = p / , -x *(Pio 9) s(p 7 , p c ) [s(g, p e ) t(p e , fc) + s(g, p 7 ) *(p 7 , fc)] (A2) 

+ q) 



M b LL 



M b LR 



M C LL 



M LR 



M d LL 



M d LR 



1 



Pw(Pe~Pu) Pw(q + q) 

x{ -2f(p„,g) s(q,p e ) [*(p 7 ,p e ) s(p e ,k) -f(p 7 ,p„) s(p u ,k)} 

+ (1 + k + X(p e -Pu) 2 )t(q,p,,) s(k, q) t(p v ,pj s(p 7 ,p e ) 

-(1 + k + \(q + q) 2 )t(p„,Pj) s(fc,p e ) *(g,P 7 ) S (P7>9) 

+\[t(p 1 ,p e ) s(p e ,k) -f(p 7 ,p„) s(p„,A;)] 

xt(g,p 7 ) s(p 7 ,g) t(p„,p 7 ) s(p 7 ,Pe)} (A3) 
1 

Pw(p e -Pu) Pw(q + q) 

x { - 2t(p„, g) s(g, p e ) [s(p 7 , p e ) * (p e , fc) - s(p 7 , p„) t(p„, fc)] 
+ (1 + re + A(p e -Pu) 2 )t(q, fc) s(p 7 ,g) t(p„,p 7 ) s(p 7 ,Pe) 
-(1 + k + A(g + g) 2 )f (p„, A;) s(p 7 ,p e ) t(g,p 7 ) s(p 7 ,g) 

+ A[s(p 7 ,p e ) *(p e , fc) - S(p 1 ,p v ) t(p„, fc)] 

xt(g,p 7 ) s(p 7 ,g) t(p„,p 7 ) s(p 7 ,p e )} (A4) 



1 (p 7 


- <?) 2 iV(Pe 


-Pv) 










Xt(q,Pu) 


s(k, q)[s(p e ,q) t(g,p 7 ) + 


s(Pe,Pv) t(p v ,Pi)] 


(A5) 




2Q f 










1 (p 7 


- q) 2 P w (p e 


~Pv) 










xt(q,Pu) 


s(p 7 ,g)[ s (Pe ? <?) t(q, k) + 


s(Pe,Pv) t(p„,k)] 


(A6) 




2Q f 










V 


-p 7 ) 2 iV(Pe 


~Pv) 










xt(g,p 7 ) 


s(p e ,q)[s(k,q) t(q,p v ) - 


s(fc,P 7 ) *(p 7 » 


Pv)] 


(A7) 


+ 7- 


._ \9 r> 


.t(q, k) s(p e ,q) s(p 7 


,q) t(q,p v ) 




(A8) 



where the propagators are defined by 

Pwip) = [P 2 ~M^ + iT w M w }. (A9) 
The first subscript of the amplitudes refers to the helicity of the electron and the second 
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subscript to the helicity of the photon. The amplitudes correspond to the diagrams of 
Fig. 2 where the four momenta p e ,p 1 ,p u ,q and q are defined. Qf = Qd = ~ 1/3 and 
Qf = Qu = +2/3. Helicity amplitudes not explicitly written down are zero. To obtain the 
cross section the amplitudes for given electron and photon helicities are summed over and 
squared. These are then averaged to obtain the spin averaged matrix element squared and 
finally integrated over the final state phase space to yield the cross section. 

APPENDIX B: THE PHOTON DISTRIBUTIONS 

To obtain the cross sections in the main text we convolute the ej cross section with the 
relevant photon distributions: 

cr = J fj/ e (x) a(e7 — > Wv) dx (Bl) 

where the various photon distributions, / 7 / e (a;), are given below. 



1. Back-Scattered Laser Photons 



Intense high energy photon beams can be obtained by backscattering a low energy laser 
off of a high energy electron beam. The energy spectrum of the back-scattered laser photons 
is given by |21 



1 



D(0 



1 - x + 



1 



Ax Ax 2 
+ 



(B2) 



l-x f(l-z) i 2 {l-x) 2 _ 

where the fraction x represents the ratio of the scattered photon energy, u, and the initial 
electron energy, E, (x = uj/E) and 

1 



D(0 



with 



\ 4 8 \ , /-, ^ 1 8 



4Eu 2 a 2y/suj 
f = — cos — ^ — 

?7ig 2 m A e 



(B3) 



(B4) 
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and ljo is the laser photon energy and «o ~ is the electron-laser collision angle. The 
maximum value of x is 

Xm = IT = (TTO ■ (B5) 

Because of the onset of e + e~ pair production between backscattered and laser photons, 
conversion efficiency drops considerably for x > 2 + 2y2 ~ 4.82. We use this value which 
for 250 GeV electrons corresponds to a laser energy of about 1.26 eV. 

The photon spectrum is sensitive to the the product A e P 7 where A e is the mean elec- 
tron helicity and P 7 is the mean laser helicity. Larger values of A e P 7 give a harder more 
monochromatic photon spectrum. Measuring the actual A e introduces systematic errors so 
we assume that the electron beam is unpolarized. In constrast the laser can be easily polar- 
ized almost completely. The amount of polarization is energy dependent. Assuming P 7 = 1 
the average helicity £2 of the photon beam is given by 

, = at-2x-tx)(2-2x + x 2 ) 

^ 2 2£ 2 - 4x£ - 4£ 2 x + Ax 2 + 4£x 2 + 3£ 2 x 2 - £ 2 a; 3 1 ' 

The long dashed line in Fig. 13 shows the spectrum of photons for an unpolarized laser and 
the medium dashed line shows the photon spectrum with helicity — P 7 . Note that the most 
energetic photons are always polarized with opposite helicity to that of the laser photons. 

2. Beamstrahlung Photons 

The interpenetration of the dense electron and positron bunches in future e + e~ colliders 
generates strong accelerations on the electrons and positrons near the interaction point which 
gives rise to synchrotron radiation which is referred to in the literature as beamstrahlung 
|23|p4]] . Beamstrahlung depends strongly on machine parameters such as luminosity, 



pulse rate, and bunch geometry. The distribution function of beamstrahlung photons can 
be written in the following approximation: 

ff/T{x, b) = f^ e (x, b) &(x c -x) + f$(x, b) Q(x - x c ) (B7) 
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where, as before, x is the fraction of the beam energy carried by the photon, b is the impact 
parameter of the produced 7, and x c separates low and high photon energy regions where 
the different approximations to f^ m (x,b) are used. The distribution used for small and 
intermediate values of x is given by 



4/1 M) 



CK 
TVs 



(1-x) 



X 



2/3 d 



X 



1 1/3 



1 + 



6CT 2 / 3 VI 



,r 



2/3 



exp 



x 



3T 



:bs) 



where C = —Ai'(0) = 0.2588, and Ai(x) is Airy's function. For large values of x the 
distribution is given by 



A" 



2 v ^FT 1 /2 



rrfl — a;) 



;l/2(l_a;)l/2 



exp 



3T 1-x 



(B9) 



The value x c is such that /,^ m is continuous at x = x c . It depends on the machine design; 
for the 500 GeV NLC, x c ~ 0.48. The dimensionless quantities A and T are defined as 

a z E ± 



A = 2y/3. 



a- 



m 



T 



(BIO) 



where m and p are the electron mass and momentum, and E± is the transverse electric field 
inside a uniform elliptical bunch of dimensions l x y = 2cr x ,y an d l z = 2^/3a z 



Na 



1,2 \ 1/2 



+ 



(Bll) 



with A being the number of particles in the bunch. 

For the case of beamstralung we have to average over the impact parameter in addition 
to integrating over the energy fraction, x. Using elliptical coordinates the expression in the 
brackets of Eq. pil| becomes 



L2 \ 1/2 

y 



(B12) 



and 



fin?) 



f f b jr(x,b)2bdb 



(B13) 
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The photon luminosity of beamstrahlung is very sensitive to the transverse shape of the 
beam. The aspect ratio 

G = (B14) 

provides a good measure of beamstrahlung with large photon luminosities associated with 
small values of G. For high photon luminosity one tunes to round beams, G — 1. For the 
original NLC design G ~ 2.7. We include the beam parameters for a number of NLC options 
in Table III. 




x ( = w/E) 

Figure 13: The photon distributions described in the text. The solid line is the Weizacker- 
Williams distribution, the long dashed line is the backscattered photon distribution, the medium 
dashed line is the backscattered photon distribution for photons with opposite polarization to that 
of the laser, the short dashed line is the backscattered photon distribution for photons with the 
same polarization as the laser, the dotted line is the beamstrahlung spectrum for y/s = 500 GeV 
with G=2.7, and the dash-dot line is with G=l. 
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TABLE HI. NLC Machine Parameters 

NLC TLC 

E cm (TeV) 0.5 0.5 1 1 
L (cm^sec- 1 ) 9 x 10 33 

N 1.67 x 10 10 1.67 x 10 10 1.67 x 10 10 1.67 x 10 10 

a z (cm) 0.011 0.011 0.011 0.011 

o y (cm) 6.5 x 10~ 7 3.3 x 10~ 6 1.7 x 10" 5 3.3 x 10" 6 

a x (cm) 1.7 x 10~ 5 3.3 x 10~ 6 1.7 x 10" 5 3.3 x 10" 6 

G 2.7 1.0 2.7 1.0 

x r 0.48 0.64 0.58 0.72 



3. Classical Bremstrahlung 

Finally we consider conventional bremsstrahlung of photons by electrons which we use 
for the a/s = 200 GeV case and which also contributes to the photon luminosity when 



we consider beamstrahlung. We use the well-known Weiszacker- Williams distribution f44 
which we include for completeness: 



fiSTix, E„) = £ 1 + ( V" )J In (%- | (B15) 
where E max is the electron beam energy. 
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